Phase diagram of the t — t' Hubbard model taking into account spin-spiral waves and 

phase separation at finite temperatures 
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The effect of temperature on the magnetic phase separation and the parameters of spin-spiral 
waves (SSW) is studied using a two-dimensional (2D) single-band t — t' Hubbard model and the 
Hubbard- Stratonovich transformation. Both commensurate (antiferromagnetic (AF)) and incom- 
mensurate (helical) magnetic phases are considered. It is shown that the temperature significantly 
affects the collinear and helical magnetic phases. With an increase in the temperature, the phase- 
separation (PS) regions (AF+ [Q,Q]), ([Q,Q] + [Q,t]) get substantially reduced but new regions 
([Qi,tt] + [Q2,7r]), (AF+[Q, 7t]) arise. The results are used for the interpretation of the magnetic 
properties of cuprates. 
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The discovery, in the late 1980s, of high-temperature 
superconductivity giving rise to new extensive studies of 
the phenomenon of superconductivity, has also regener- 
ated a vivid interest in quantum low-dimensional mag- 
netism. At present most researchers are inclined to be- 
lieve that magnetic interactions are responsible in great 
part (if not completely) for the high-temperature super- 
conductivity. Apparently, the high-temperature pseudo- 
gap observed in superconductors is also connected with 
magnetism, rather than with the Cooper pairs [1]. The- 
oretical description of magnetism in these systems is still 
far from complete. As is shown experimentally [2], the 
itinerant magnetism should enter as an essential compo- 
nent into theoretical models. Of great importance there- 
with is the two-dimensionality [3, 4]. The authors of 
paper [5] have studied the ground-state phase diagram 
of the 2D single-band Hubbard model. The calcula- 
tions were performed in the framework of the mean field 
approximation with consideration of both commensu- 
rate (ferromagnetic and AF) and incommensurate (SSW) 
magnetic states. The results obtained made it possi- 
ble to explain some peculiarities of the behavior of su- 
perconducting cuprates, and clarified the role of spatial 
magnetic phase separation (PS). Further investigation of 
the model at finite temperatures [6] has shown the ne- 
cessity of taking into account the next-nearest-neighbor 
hopping, and revealed a significant influence of thermo- 
dynamic fluctuations on the boundaries of PS regions. 

In this paper we consider the 2D single-band Hubbard 
model at finite temperatures with allowance for the next- 
nearest-neighbor hopping. Note a controversial point 
that arises in studying the 2D models. It is related to 
the well known statement that magnetic order is lack- 
ing in 2D systems at finite temperatures (except for the 
Ising-type models). As a rule, the researchers which use 
the 2D Hubbard model for describing cuprates assume 
that the weak magnetic coupling between the Cu lay- 
ers suppresses anomalous growth of transverse magnetic 
fluctuations and stabilizes the order, ensuring a quasi- 



two-dimensional behavior of the magnetic characteristics. 
However there may be an alternative mechanism of sup- 
pressing the anomalous growth of magnetic fluctuations 
through the charge degree of freedom and the volume 
Coulomb interaction. In this case one can speak of a 2D 
magnetic subsystem. The problem of the existence of 2D 
magnetism is similar to that of the 2D crystals in which 
the anomalous growth of lattice fluctuations is suppressed 
owing to strong anharmonicity of the longitudinal and 
transverse lattice vibrations [7]. In view of the above, in 
this work we use a special approximation which makes it 
possible to prevent anomalous growth of magnetic fluctu- 
ations and stabilize the magnetic order of the 2D model, 
leaving the question of the existence of 2D ordered struc- 
tures unsettled. It should be noted that the results of 
such approximations are useful even in the absence of 
long-range magnetic order, as they may be used for the 
description of the dominant thermodynamics! magnetic 
fluctuations at fixed temperatures. Another controversial 
point of our study is the choice of possible incommensu- 
rate magnetic structures. We have restricted our con- 
sideration to the helical magnetic phases (SSW), leaving 
aside the collinear spin-density waves (SDW). Usually, 
the findings obtained in neutron experiments give no pos- 
sibility to distinguish between SSW and SDW (see, e.g., 
[8-11], except in a few cases such as the spin-polarized ex- 
periments on Fe-As monocrystals [12] where the existence 
of SSW has been established unambiguously. Ab initio 
calculations show that SSW may be realized in many 
systems on basis of the transition metals [13]. A direct 
comparison of the SSW and SDW energies within the 2D 
Hubbard model shows that in some cases the SDW are 
energetically preferable [14]. At the same time the ob- 
tained charge redistribution over the sites in SDW is too 
large, and the Coulomb interaction, not accounted for 
in the model, should increase the SDW energy, making 
these states unfavorable. 



2 



MODEL 

We adopt the Hamiltonian of the Hubbard model on 
the squared lattice H = Hq + Hi nt , 

Ho = t^tsCjsi Hint = U^2 %t%> (1) 

ijs j 

where tij = —t for the nearest-neighbor sites i,j, tij = t' 
for the next-nearest neighbors, hj s = c^ s Cj s is the elec- 
trons number operator with spin projection s =t,4- on 
site j, c~j s (Cj s ) is a creation (annihilation) electronic op- 
erator with spin projection s on site j. U denotes the 
interaatomic Coulomb interaction on site. 

We consider the helical-type magnetic structures cor- 
responding to the magnetization vector uniformaly ro- 
tating in the polarization plane when moving from one 
site to another. This magnetic structure is character- 
ized by the magnitude and direction of the wave vector 
Q = [Qx , Qy] ■ Generally, the wave vector Q does not co- 
incide with a reciprocal lattice vector and turns out to be 
incommensurate [5, 6]. The superposition of the helical 
wave and the ferromagnetic component perpendicular to 
the polarization plane has a large energy [15] and it is 
not considering here. 

The thermal properties of the system are determined 
by the partition function of the grand canonical ensemble 









Z = Tr 


T T cxp 


-jT P H(r)dr| 









Here the symbol Tr denotes summation over the comlete 
set of quantum states, T T denotes the time ordering oper- 
ator, (3 = 1/ksT denotes the inverse temperature, H(t) 
is the H — /in operator in the interaction representation 
where (i and n are the chemical potential and the total 
electron number operator, respectively. 

Partition function (2) for the many-particle problem 
is transformed into that of the single-particle problem 
with time-dependent fictitious fields, by means of the 
functional integral method through the use Hubbard- 
Stratonovich transformation [16]. To reproduce the re- 
sults of the generalized Hartree-Fock approximation at 
the ground state [5, 6] we need to use the two-field func- 
tional integral method in the static approximation [17] 
introducing the spin and charge auxiliary fictitious fields, 
Vj and Q. In addition, we adopted the saddle point 
approximation for the charge auxiliary fictitious field Q 
what corresponds to the replacement of the Q by its value 
C°(vj) at the saddle point. Using C^( v j); we minimize the 
thermodynamic potential at the fixed configuration Vj, 
but neglect the charge fluctuations. In this work, we take 
into account only the longitudinal spin fluctuations spec- 
ifying the direction of the spin auxiliary fictitious field Vj 
parallel to the magnetization vector m., at each site. 

The partition function is expressed by the mean field 
approximation for the single-particle Green function 



G(z) = [z — Mmf] -1 [18, 19], where Hmf is the mean 
field approximation for the Hubbard Hamiltonian (1). 
Finally, we introduce the self-energy of electron in the 
effective medium S, which is determined by the average 
with respect to the amplitude of spin fictitious field v 
Green function G =< [1 - GViv^G >. The effective 
self-energy £ and G are found using the self-consistent 
matrix equation in the single-site coherent potential ap- 
proximation 

G(z)=<[l-GV(v)]~ 1 G >=G(z-E). (3) 

After introducing these apprpximations, the partition 
function is represented in the form of an integral over 
the amplitude of spin fictitious field v 

Z = cxp{-/3(ft[£] + AO)}, (4) 
fi[£] = VLmf — -q lndet[l + £(ioj„)G(iu; n )], 

n 

An = -hn { dvexp{~/3An(v)}, 

P J 

Afi(u) = ^-^2lndet[l - G(iu n ){V{v) - E(iw„))] 

n 

+U [v 2 + ( 2 (v)] . 

Here fi[£] is the thermodynamic potential for the 
effective medium, which is determined by the ef- 
fective self-energy S, and Af2 is the fluctuat- 
ing part of the thermodynamic potential, Ojv/f = 
— l//31nTrexp [— /3(Hmf ~ A*^)] is th e mean field ap- 
proximation for the thermodinamical potential of elec- 
trons, u> n = 7r(2n + l)//3 are the Matsubara frequencies 
for the Fermi particles. For brevity, site subscript in (4) 
is omitted. 

We mentioned earlier that the ground state in this ap- 
proximation is the Hartree-Fock state, which is formally 
validated in the limit of the small parameter U/t. In our 
theoretical approach, the analysis of the approximations 
in terms of the small parameter at nonzero temperature 
is an extremely complicated problem, which is still un- 
solved. Nevertheless, the static approximation is justified 
since we are not studying the superconducting properties 
of our systems and the strong correlations in the chosen 
range of parametrs at the electron density close to unity 
are not very significant. This confirmed by comparison 
of our results [5] with those obtaned by the dynamical 
claster methods [20]. 

In spite of the aforementioned essential approxima- 
tions, further analysis of expression (4) is possible only by 
using numerical methods. The self-consistent solution to 
(3) and (4) allows us to calculate all magnetic properties 
of our system under the condition that the magnetic state 
with the minimum thermodynamic potential is chosen. 
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RESULTS 

All calculations were carried out at U/t = 4.8, t' /t = 
0.2. For these values various magnetic phases are re- 
alized [5], and the characteristics of the electron states 
of superconducting cuprates (La,(2- x )Sr x Cu04) [21] may 
be described well enough. We restricted ourselves to the 
electron concentrations less than unity, because in a wide 
range of larger concentrations the system exhibits a sta- 
ble AF ordering [">]. Figure 1 presents the T — n phase 
diagram (hereafter we shall use a dimensionless temper- 
ature T — > ksT/t). Thick solid lines correspond to the 
phase transitions between states with different magnetic 
order. The phase transition from the ordered magnetic 
state to a paramagnet (P) is a second-order transition, 
all the other are first-order transitions. This is well il- 
lustrated by the dependence of chemical potential on the 
electron concentration at T = 0.02 Fig. 2 which exhibits 
a inflection when passing from AF to the diagonal [Q, Q] 
phase, and from [Q, Q] to the parallel [Q, ir] phase. In- 
stability of the chemical potential (negative derivative 
with respect to concentration) results in spatial magnetic 
phase separation (PS) whose boundaries are determined 
by the Maxwell rule: 

/ Km) - M(n)] dn = 0, (5) 

Jnl 

where nl and n2 are the boundaries of the PS region. A 
comparison of the PS region boundaries at zero and finite 
temperatures shows that with increasing temperature 
the PS regions narrow down, being replaced by homo- 
geneous states. Besides there arise three new PS regions: 
{[QiM + [Qs.tt]), ([Q,Q] + [Q,tt]), and (AF+[Q,n]) 
which were lacking at zero temperature. The charac- 
teristics of the magnetic phases and their partial ratios 
inside the PS regions are completely determined by their 
values at the boundaries of these regions. Figures 3, 4, 5, 
6 present the parameters of spiral magnetic structures 
at the PS boundaries (AF+ [Q,Q]), ([Q, Q] + [Q, it]), 
([<9i,7t] + [Q 2 ,tt]), and (AF + [Q,7rj). These plots may 
be used to determine the mean local magnetization at 
site < ?7i > , the mean absolute magnitude of local mag- 
netization < \m\ >, the spiral wave vector Q, and the 
partial phase ratio for any point inside the PS region. 

Of interest is a rather strong T dependence of Q in 
the phase mixture (AF+ [Q,Q]), see Fig. 3. In experi- 
ments on cuprates close to half filling a region, generally 
referred to as "spin glass" (SG), is observed [4]. The PS 
(AF+ [Q,Q]) with a strong temperature dependence of 
vector Q found in our study can be easily associated with 
the experimentally observed SG region. Note that the ra- 
tio between the temperatures T/y and T g corresponding 
to the P— > AF and AF— > SG transitions, respectively, 
for n = 0.98 is approximately equal to 12 [11], which 
agrees with the ratio of temperatures T/v and Tps for 



the P — > AF and AF — > (AF + [Q, Q}) transitions, respec- 
tively, in our calculation Fig. 1. Recall that the calcu- 
lated absolute values of the transition temperatures are 
overestimated because of the model and computational 
approximations made (see Section 1). The distinctive 
feature of the PS regions ([Qi,7r] + [Q2,tt]) Fig. 5 and 
(AF + [Q, 7r]) Fig. 6 consists in a large difference between 
the mean local magnetic moments in the phase mixture, 
the absolute magnitudes of local magnetic moments be- 
ing approximately equal. At high temperatures this fact 
can be interpreted as the effect of long-lived fluctuations 
against a paramagnetic-state background [3], Besides, 
the superparamagnetic behavior experimentally observed 
in chemically homogeneous Fe-Al alloys at high tempera- 
tures [22] is easily explained when taking into account the 
possibility of spatial magnetic separation between SSW 
and the ferrimagnetic phase. Note that low-temperature 
neutron experiments have revealed in these alloys the ex- 
istence of SDW [23] the parameters of which are fairly 
well described in ab initio calculations using SSW [22, 24]. 

The temperature and concentrational behavior of the 
PS region ([Q,Q] + [Q,7r]) reproduces two significant 
experimental facts [10]. A decrease in the number of 
electrons results in the following sequence of transi- 
tions: [Q,Q] -> {[Q,Q] + [Q,tt]) -> [Q,tt], and an 
increase in the temperature gives rise to the transition 
([Q,Q] + [Q,tt]) -)• [<3,tt], see Figs. 1, 4. It should be 
noted that quantitative agreement with the experimen- 
tal concentrations of the transitions cannot be expected 
in our calculations, as they are extremely sensitive to the 
parameters U/t and t' /t. An insignificant decrease in U/t 
from 4.8 leads to a considerable shift of these transitions 
towards concentrations close to unity, and at U/t < 4 
the concentrational transition connected with the [Q, Q] 
phase totally disappears, see [5]. 

The authors of paper [5] have pointed out the pro- 
portional dependence of the wave vector on the electron 
concentration, which is in agreement with the experiment 
[10] and can be easily interpreted in terms of the spin- 
spiral structure formation. At finite temperatures this 
dependence is retained, see, e.g. Fig. 2. Moreover, in the 
considered range of concentrations n, a general tendency 
for an increase in the SSW wave vector with tempera- 
ture is observed. This is true both for each of the PS 
phases Figs. 3-6 and for the homogeneous states Fig. 7. 
However note that in the PS region the mean wave vec- 
tor can diminish with increasing temperature. This is 
the case, for example, in the PS ([<3i,7r] + [<52,tt]) be- 
cause of a change in the partial phase ratio, namely, a 
decrease with increasing temperature in the portion of 
the phase with larger wave vector Fig. 5. Of particular 
interest is the temperature transition metal-isolator near 
the concentration equal to unity. 

Figures 8 and 9 show the density of states for the AF 
and [Q, Q] phases in the PS (AF+ [Q, Q]) region at tem- 
peratures 0.02 and 0.002. It is seen that for T = 0.02 all 
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boundaries of PS AF+(Q,Q) 
(Q,Q)+(Q,Jt) 
(Q„7t)+(Qj,Jt) 
AF+(Q,7C) 
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Figure 1: (Color online) Magnetic phase diagram at U/t = 4.8 
and t! /t = 0.2. Blue bold line denote the Neel tempera- 
ture (second-order phase transition), black bold lines denote 
first-order phase transitions calculated without regard for PS, 
shaded areas denote PS regions. 



Figure 2: (Color online) Chemical potential fi (left axis) and 
wave vector Q (right axis) versus the electron density n at 
U/t — 4.8 and t' /t = 0.2. Dashed lines denote n dependence 
of the [i in the PS regions. 



the phases are metallic, and AF phase may be considered 
as a conductor, since it has a lower electron density at 
the Fermi level as compared to the [Q, Q] phase Fig. 8(a). 
With decreasing temperature the AF phase actually be- 
comes isolator, and the conductivity of the [Q, Q] phase 
increases Fig. 9(a). The metal-isolator transition experi- 
enced by the AF phase is a Slater-type transition [25]. 
The total density of states for temperatures 0.02 and 
0.002, and concentration n = 0.96 is shown in Fig. 8(b) 
and Fig. 9(b). In both cases the total density of elec- 
tron states at the Fermi level is different from zero. It is 
obvious however that at low temperature the system is 
practically an isolator, because at this temperature and 
concentration the portion of the AF phase relative to the 
phase exceeds 1/2 (the percolation threshold in 
the 2D case being equal to 1/2). Thus our system will 
undergo a metal-isolator transition in the concentration 
range 1 0.95, but this is a Slater-type transition of per- 
colation character. We did not account for the strong 
electron correlations which could modify the boundaries 
and features of the region of the transition described. It 
is also evident that with a decrease in the electron con- 
centration (hole doping) the role of correlations substan- 
tially grows. The Mott transitions induced by the strong 
correlations in optimally doped and overdoped cuprate 
regions have been considered in detail in the review [26]. 



CONCLUSION 

Thus, based on 2D Hubbard model, we have studied 
the formation of helical waves (SSW) and the magnetic 
PS at finite temperatures. The calculated T — n phase 
diagram Fig. 1 shows that the PS regions existing at 
T = with increasing temperature get much narrower, 
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Figure 3: (Color online) Electron number dependence of the 
average local magnetic moment < m > and its absolute value 
< \m\ > (left axis), wave vector Q (right axis) along the 
boundary of AF + [Q, Q] PS region at U/t = 4.8 and t' /t = 
0.2. Dashed lines - [Q, Q] phase, bold lines - AF. 
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Figure 4: (Color online) [Q,Q] + [Q,7r] PS region. Dashed 
lines - [Q,Q] phase, bold lines - [Q,tt]. Notations are the 
same as in Fig. 3. 



and there arise new PS areas with different symmetry of 
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Figure 5: (Color online) [Qi,7r] + [Qzi""] PS region. Dashed 
lines - top boundary whith [Q, n] phase, bold lines - lower 
boundary whith [Q,ir] phase. Notations are the same as in 
Fig. 3. 
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Figure 8: (Color online) Electron density of states EDS in the 
AF + [Q, Q] PS region, (a) Porcion EDS on the boundaries 
of PS at T = 0.02. (b) Total EDS at T = 0.02 and n = 0.96. 
Dashed line - Fermy level U/t = 4.8, t! /t = 0.2 
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Figure 6: (Color online) AF+ [Q,n] PS region. Dashed lines 
- AF, bold lines - [Q, n]. Notations are the same as in Fig. 3. 
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Figure 7: (Color online) Temperature dependence of the wave 
vector Q (left axis), average local magnetic moment < m > 
and its magnitude < |m| > (right axis) at U/t = 4.8, t' /t = 
0.2 and n = 0.88. 



the wave vectors SSW. The behavior of the wave vec- 
tors, mean local magnetic moments, and mean absolute 
magnitudes of the local moments is described for both 



the homogeneous SSW states and the phases forming 
the PS region. It is shown that the mean local mag- 
netic moments of different phases inside the PS region 
may considerably differ in magnitude and temperature 
behavior. The wave vectors for all SSW increase with 
temperature. The proportionality of the SSW wave vec- 
tors to the number of electrons (holes) is retained in a 
wide range of temperatures. 

Taking account of the model and mathematical ap- 
proximations made, the results obtained may be used 
to describe quasi-two-dimensional systems or the param- 
eters of thermal magnetic fluctuations in the range of 
electron concentrations and temperatures where the ef- 
fect of the strong correlations and system dynamics is 
not too significant. 

A comparison of the results obtained with the avail- 
able experimental data for high-temperature compounds 
La(2-x)S r x CuOi shows a good semiquantitative agree- 
ment as to the behavior of the magnetic characteris- 
tics in the half-filling region with small hole concentra- 
tion. The ratio between the temperatures (P — > AF 
transition) and T g (AF — > "spin glass" transition) at 
n = 0.98 approximates to 12 [11], which agrees with 
the ratio of temperatures Tn and Tps for the P — > 
AF and AF— > (AF+ [Q, Q]) transition, respectively, ob- 
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Figure 9: (Color online) (a) T = 0.002, n = 0.9044 and n = 
0.9970. (b) T = 0.002 and n = 0.96. Notations are the same 
as in Fig. 8. 



tained in our calculation. The presence of the PS region 
([QjQ] + [Q,tt])) the sequence of transitions [Q,Q] — > 
([Q, Q] + [Q, tt]) — > [Q,tt] occurring with increasing con- 
centration, and ([Q, Q] + [Q,tt]) — >• [Q,7r] taking place 
with increasing temperature coincide with the PS region 
([Qj<9] + [Q,7r]) and the transition sequences observed 
experimentally [10]. 

Taking into account the possibility of a great difference 
in the mean local magnetic moment between the phases 
forming the PS region, one can explain the superparam- 
agnetic behavior of chemically homogeneous alloys Fe-Al 
[22]. 

Obviously, allowance for the strong correlations and 
the system dynamics will affect the quantitative char- 
acteristics of the PS regions and the SSW parameters. 
However, as demonstrated by a comparison of TV ob- 
tained in our calculations and in those performed within 
the dynamic cluster approximation (DC A) [20], the dif- 
ference proves to be not very large. 

This work was supported by the Russian Foundation 
for Basic Research (project No. 12-02-00632) and Ur Br 
RAS (project No. 12-U-2-1021). 



* Electronic address: arzhnikof@bk.ru 
' Electronic address: groshev_a.g@mail.ru 
[1] Rui-Hua He, M. Hashimoto, H. Karapetyan, J.D. Ko- 
ralek, J. P. Hinton, J. P. Testaud, V. Nathan, Y. Yoshida, 
Hong Yao, K. Tanaka, W. Meevasana, R.G. Moore, 

D. H. Lu, S.-K. Mo, M. Ishikado, H. Eisaki, Z. Hussain, 
T.P. Devereaux, S.A. Kivelson, J. Orenstein, A. Kapit- 
ulnik, Z.-X. Shen, Science 331, 1579 (2011). 

[2] Rui-Hua He, M. Fujita, M. Enoki, M. Hashimoto, 

S. Iikubo, S.-K. Mo, Hong Yao, T. Adachi, Y. Koike, 

Z. Hussain, Z.-X. Shen, and K. Yamada, Phys. Rev. Lett. 

107, 127002 (2011). 
[3] Shirane, G.Y. Endoh, R.J. Birgeneau, M.A. Kast- 

ner,Y Hidaka,M. Oda,M. Suzuki, and T. Murakami, Phys. 

Rev. Lett. 59, 1613 (1987). 
[4] M.A. Kastner, R.J. Birgeneau, G. Shirane ,Y. Endoh, 

Rev. Mod. Phys. 70, 897 (1998). 
[5] P.A. Igoshev, M.A. Timirgazin, A.A. Katanin A.K. Arzh- 

nikov, and VYu. Irkhinl, Phys. Rev. B 81, 094407 

(2010) . 

[6] A.K. Arzhnikov, A.G. Groshev, JETP Lett. 94, 703 

(2011) . 

[7] A. Fasolino, J.H. Los and M.I. Katsnelson, Nature 78, 
2011 (2007). 

[8] K. Yamada, C.H. Lee, K. Kurahashi, J. Wada, S. Waki- 
moto, S. Ueki, H. Kimura, Y. Endoh, S. Hosoya, G. Shi- 
rane, R.J. Birgeneau, M. Greven, M.A. Kastner, and 
Y.J. Kim, Phys. Rev. B 57, 6165 (1998). 

[9] M. Matsuda, M. Fujita, K. Yamada, R.J. Birgeneau, 
M.A. Kastner, H. Hiraka, Y. Endoh, S. Wakimoto and 
G. Shirane, Phys. Rev. B 62, 9148 (2000). 
[10] M. Fujita, K. Yamada, H. Hiraka, P.M. Gehring, 
S.H. Lee, S. Wakimoto and G. Shirane, Phys. Rev. B 65, 
064505 (2002). 

[11] M. Matsuda, M. Fujita, K. Yamada, R.J. Birgeneau, 
Y. Endoh and G. Shirane, Phys. Rev. B 65, 134515 
(2002). 

[12] E.E. Rodriguez, C. Stock, K.L. Krycka, C.F. Majkrzak, 
P. Zajdel, K. Kirshenbaum, N.P. Butch, SR. Saha, 
J. Paglione, and M.A. Green, Phys. Rev. B 83, 134438 
(2011). 

[13] R. Lizarraga, L. Nordstrom, L. Bergqvist, A. Bergman, 

E. Sjostedt, P. Mohn and O. Eriksson, Phys. Rev. Lett. 
93, 107205 (2004). 

[14] M.A. Timirgazin, A.K. Arzhnikov and A.V Vedyaev, 
Solid State Phenom. 190, 67 (2012). 

[15] M.A. Timirgazin and A.K. Arzhnikov, Solid State Phe- 
nom. 559, 152 (2009). 

[16] Yu.A. Izyumov and Yu.N. Skryabin, Statistical Mechan- 
ics of Magnetically Ordered Systems (Nauka, Moscow, 
1987; Springer- Verlag, New York, 1988). 

[17] R.F. Hassing, D.M. Esterling, Phys. Rev. B 7, 432 
(1973). 

[18] S.Q. Wang, W.E. Evenson, and JR. Schrieffer, Phys. 

Rev. Lett. 23, 92 (1969). 
[19] A.A. Abrikosov, LP. Gor'kov, and I.E. Dzyaloshinskii, 

Methods of Quantum Field Theory in Statistical Physics 

(Fizmatlit, Moscow, 1963; Dover, New York, 1975). 
[20] T. Maier, M. Jarrell, T. Pruschke, M.H. Hettler, Rev. 

Mod. Phys, 77, 1027, (2005). 
[21] A.I. Lichtenstein, M.I. Katsnelson, Phys. Rev. B 62, 



7 



R9283 (2000). 

[22] A.K. Arzhnikov, L.V. Dobysheva, M.A. Timirgazin, J. 

Magn. Magn. Mater. 320, 1904 (2008). 
[23] D.R. Noakes, A.S. Arrott, M.G. Belk, S.C. Deevi, 

Q.Z. Huang, J.W. Lynn, R.D. Shull, D. Wu, Phys. Rev. 

Lett. 91, 217201 (2003). 
[24] J. Bogner, W. Steiner, M. Reissner, P. Mohn, P. Blaha, 



K. Schwarz, R .Krachler, H. Ipser, B. Sepiol, Phys. Rev. 

B 58, 14922 (1998). 
[25] M. Imada, A. Fujimori and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998). 
[26] Patrick A. Lee, Naoto Nagaosa and Xiao-Gang Wen, Rev. 

Mod. Phys. 78, 17 (2006). 



